Confinement-driven translocation of a flexible polymer 
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We consider the escape of a flexible, self-avoiding polymer chain out of a confined geometry. 
By means of simulations, we demonstrate that the translocation time can be described by a sim- 
ple scaling law that exhibits a nonlinear dependence on the degree of polymerization and that is 
sensitive to the nature of the confining geometry. These results contradict earlier predictions but 
are in agreement with recently confirmed geometry-dependent expressions for the free energy of 
confinement. 
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Translocation through a nanopore is one of the fun- 
damental biological mechanisms through which long 
molecules can be exchanged between different regions 
compartmentalized by biological membranes 0. Ex- 
amples of this phenomenon include the injection into 
host cells of DNA packed inside virus capsids 0] and 
the transport of proteins through biological membranes. 
Furthermore, pioneering experiments have demonstrated 
that DNA can be translocated through a nanopore by 
means of an external electric field, and that this event 
can be probed by measuring the variation in ionic cur- 
rent through the pore 0, 5- This has opened the 
prospect of creating efficient and economical DNA se- 
quencing devices and has resulted in a widespread the- 
oretical and experimental interest in polymer transloca- 
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The passage of a flexible chain through a narrow open- 
ing involves a large entropic barrier, so that most poly- 
mer translocation phenomena require a driving force. A 
typical experimental setup |2| consists of two chambers 
separated by an interface. A DNA molecule is placed in 
one of the chambers and passes through the interface via 
a small orifice, either a protein complex embedded in a 
membrane or a solid state nanopore. In this case, the 
required driving force is provided by an external electric 
field. However, one can also envisage the use of other 
forces, e.g., generated by optical tweezers or by an os- 
motic pressure resulting from the geometrical confine- 
ment of the polymer in one of the chambers. 

In the study of polymer translocation, the duration of 
the sequential passage of a chain through a membrane — 
measured from the entrance of the first monomer into 
the pore — occupies a central place, because it is one of 
the few dynamical para meters that is accessible to cur- 
rent experiments |3.ll5l Hfl Ht| . In particular, there have 
been various efforts to determine and to understand the 
dependence of translocation time r on the degree of poly- 
merization and the magnitude of the driving force. Sung 
and Park [j| proposed to treat polymer translocation as a 
diffusion process across a free energy barrier. Muthuku- 
mar [gj reproduced this approach using nucleation theory, 
employing a corrected diffusivity. This analysis reveals 
the existence of two distinct regimes. If the chemical po- 



tential gradient per monomer Afi is sufficiently small, the 
entropic barrier dominates the diffusion process, 



N 2 for N\A(i\ < 1 , 



(1) 



whereas stronger driving forces lead to a drift-dominated 
regime, 



N 

— for N\A(i\ > 1 
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Chuang et al. |3] observed an inconsistency in the rea- 
soning leading to Eq. J3J, since it assumes that the poly- 
mer chain is equilibrated at all times during the translo- 
cation, whereas the Rouse equilibration time scales as 
N 2lJ+1 and hence for sufficiently long chains exceeds the 
predicted translocation time. Furthermore, it was ar- 
gued [lflj that the unhindered motion of a chain pro- 
vides a lower bound for the translocation time, such that 
Eq. J2J should be replaced by 



(3) 



The difference in chemical potential of the monomers 
on either side of the interface clearly is a crucial ingre- 
dient in all estimates of the translocation time. Unlike 
other simulation studies [ill fl^| , here we concentrate on 
the situation where this difference originates from poly- 
mer confinement, i.e., where the chain is ejected out of 
a restricted geometry. This situation has been analyzed 
m Refs. @ and H, where numerical results were pre- 
sented that agree with Eq. (J2J) , rather than with the cor- 
rected prediction Eq. ©. To confuse matters further, it 
was recently established that the drivin g fo rce exhibits a 
nontrivial dependence on geometry |18 l ll9 | that was not 
taken into account properly in Refs. Q and 0- In this 
Letter we resolve this contradictory situation by means 
of an accurate numerical study of the escape of a polymer 
chain out of different confining geometries, in which we 
independently vary the degree of polymerization and the 
strength of the osmotic driving force. 

The free energy cost of confining a linear, flexible poly- 
mer within a planar or cylindrical geometry (uniaxial 
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or biaxial confinement, respectively) is given by a well- 
known blob scaling result |2Cj . 




where 7? is the separation between the plates or the ra- 
dius of the cylinder and Rq ~ crN" is the radius of gy- 
ration of the polymer in its unconfined state, a is the 
size of a monomer, v ~ 0.588 is the Flory exponent and 
/3 = l/(k&T), with fee the Boltzmann constant and T 
the absolute temperature. As has been argued on theo- 
retical grounds 0] and confirmed numerically [J^, this 
prediction is invalid for a spherical geometry (triaxial 
confinement) and must be replaced by a scaling law that 
has a much stronger dependence on the sphere radius R, 

/d \ 3/(3^-1) 
(3 AF ^(f^\ ^ ^1/(3,-1) ; (5) 

where </> = N(a/(2R)) 3 is the monomer volume fraction. 
This result can be understood within blob scaling theory 
by realizing that, unlike the planar or cylindrical case, 
the monomer concentration within a spherical cavity in- 
creases with increasing polymer size N. The extensivity 
of the free energy of confinement is then recovered only 
when a change in polymer size is accompanied by a cor- 
responding change in volume of the cavity such that the 
monomer concentration remains invariant [18l l21| . 

Since the free energy of confinement per monomer, 
AF/N, represents the chemical potential gradient that 
drives the translocation, combination of Eqs. (@J and (0 
with Eq. J2J leads to predictions for the translocation out 
of a planar or a spherical geometry, respectively, 

f R\ 1,u 

N 1+v I — ) planar confinement (6a) 

T ~ { w 

, N 1+v cf) 1/{1 - ZlJ) spherical confinement (6b) 

where the exponent 1 + v in the prefactor represents the 
lower bound proposed in Ref. To validate this pre- 
diction, we study the translocation of a flexible polymer 
chain which is modeled as a linear series of N spherical 
beads of diameter a, connected by bonds that are freely 
extensible up to a fixed value £m- All monomers interact 
via a hard-core repulsion, 

»<•*>- {I Z<°- (7) 

where is the center-to-center distance between beads 
i and j. The nearest-neighbor bonds are represented by 

Wb(r M _i) = < ' f . 8 

We mimic the dynamical properties of this model by 
means of Monte Carlo simulations in which only local, 



short-ranged displacements are employed. The monomer 
displacement per Monte Carlo step equals (Ax, Ay, Az), 
in which each cartesian component is chosen uniformly in 
the range [— 0.15tr, 0.15c] |22j. To avoid dynamical incon- 
sistencies that could result from crossing polymer bonds, 
we choose £m = V^cr. Confinement is imposed by means 
of a spherical or planar boundary of thickness <r, which 
exerts a hard-core repulsion on the monomers. For the 
planar case, the system is periodically replicated in the 
directions parallel to the plates, with a period 2Na. The 
setup is depicted schematically in Fig. ^ 




FIG. 1: Schematic setup of the Monte Carlo simulations. In 
the left-hand panel, a polymer chain is released from a spher- 
ical geometry of radius R. In the right-hand panel, the chain 
is confined between two parallel plates at separation R and 
escapes via a circular opening in one of the plates. The actual 
simulations are performed in three dimensions. 

The polymer chain is first equilibrated within the con- 
fining geometry. Subsequently, we create a smooth pore 
(shaped as the hole in a torus) of radius 1.3a to allow 
ejection of the polymer (see Fig. P). This radius is suf- 
ficiently small to practically exclude the translocation of 
folded chains. A translocation event is considered suc- 
cessful when the entire polymer escapes from the confin- 
ing geometry. We define the translocation time r as the 
difference between the time when the last monomer 
has left the cavity and the time t\ when, within a suc- 
cessful event, the first monomer has entered the pore. We 
systematically vary the degree of confinement R and de- 
termine the average escape time for chains with lengths 
ranging from iV = 40to7V = 512 monomers. For each 
choice of R and N, we generate between 400 and 1400 in- 
dependent translocation events. All simulations are per- 
formed in the drift-dominated regime, i.e., NAfi > 1, 
which requires R < Rq- 

To focus on the role of the degree of confinement, we 
calculate r for a fixed, long chain length (N = 256 for 
planar and N = 512 for spherical confinement). Figure [21 
displays r as a function of R. The escape times from 
both geometries are accurately described by a power-law 
dependence, but with strikingly different exponents. For 
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FIG. 2: Double logarithmic plot of the average translocation 
time r for a linear, flexible polymer escaping from a spherical 
and a planar geometry, as a function of the degree of con- 
finement R. The translocation times are expressed in units 
of 10 6 Monte Carlo sweeps. These data confirm the striking 
dependence on confinement geometry, Eqs. ilfinl and ijBbfl . 

confinement within a planar geometry, the driving force 
is relatively weak and we use smaller separations than for 
the spherical geometry (as small as R = 4.8a). A least- 
squares fit of both data sets to the expression t(R) = 
to + ai? 7 yields 7 = 1.54 ± 0.10 for planar confinement 
and 7 = 3.65 ± 0.08 for spherical confinement, with chi- 
square per degree of freedom (x 2 ) equal to 1.06 and 0.86, 
respectively. These results are in good agreement with 
the exponents in l/v [Eq. |@$] and 3/(3^ - 1) [Eq. 
and thus confirm the linear dependence of r on 1/A^ 
predicted by Eq. ©. The agreement is even closer if 
one notes that for the chain lengths employed here the 
effective Flory exponent is slightly larger than v = 0.588. 

Additional information about the translocation pro- 
cess can be obtained from the probability distribution 
function of translocation times P(r). Theoretical Q 
and experimental results 0, 0, El indicate that this 
distribution deviates from a Gaussian distribution and 
may be considerably skewed. Consequently, the average 
translocation time is not fully representative of the ex- 
perimental data. We sample P[f) for a chain of N = 128 
monomers, comparable to the chain length employed 
in the experimental determination of P(t) for single- 
stranded DNA . For escape from a spherical cavity, the 
driving force depends on the monomer volume fraction. 
As shown in Fig. for strong confinement (cf> = 0.15), 
the distribution is narrow and nearly Gaussian. However, 
at weaker confinement (<f> — 0.08), P(t) broadens and 
the skewness becomes clearly visible. Lubensky and Nel- 
son derived an expression that provides a reasonable 
description of the data in Fig. EJ however, this expression 
is not valid for large r. Following the experimental analy- 
sis 0| we therefore fit the data to an empirical expression 
of the form r ai exp(— a2"r). Just as in the experiments, 
the exponential term provides a good description of the 
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FIG. 3: Probability distribution P(t) of translocation 
times r, for a polymer of N = 128 monomers confined in 
a spherical cavity. The translocation times are expressed in 
units of 10 6 Monte Carlo sweeps. The narrow distribution 
corresponds to a strongly confined chain (monomer volume 
fraction <f> = 0.15) and the wider distribution is obtained 
for a smaller driving force (4> = 0.08). As discussed in the 
text, the curves represent empirical expressions of the form 
r ai exp(— <22t), with a characteristic ratio A of the width and 
the peak position that is close to typical experimental results. 

long-time tail. It is argued in Ref. that the distri- 
bution of passage times can be characterized in a useful 
way via the ratio A = <5r/r max between the width St of 
the distribution (as defined in Ref. pj) and its peak po- 
sition r max - For the distributions shown in Fig. [21 we find 
values in the range 0.53-0.56, indeed in agreement with 
experimental values A sa 0.5 0,0 and A sa 0.55 [l^ . 

We now proceed to determine the dependence of 
translocation time on the degree of polymerization N. 
For the spherical case, we perform a series of simulations 
at constant initial volume fraction <f> = 0.1, and for the 
planar geometry we perform a series of simulations at 
fixed separation R = 4.8a. For both geometries r is ac- 
curately described by a power-law dependence r ~ N s 
that is independent of geometry, in accordance with the 
observation that the free energy of confinement is 

extensive [cf. Eqs. (0J and (J5J]. For planar confinement, 
a least-squares fit yields S = 1.55 ± 0.04 (x 2 = 1.16) 
and for spherical confinement we find S = 1.59 ± 0.03 
(x 2 = 0.87). Both results agree with 1 + v i=a 1.59, con- 
firming the lower bound established by Eq. 0. Accord- 
ingly, all results for a given geometry can be combined 
in a single data collapse. Figure 0fa) (plates) shows the 
translocation time r normalized by N l+V as a function 
of the inverse driving force (A/x) -1 ~ (R/a) 1 ^, for five 
different chain lengths (N = 32, 64, 96, 128, and 256). 
Likewise, Fig. 0{b) (sphere) displays t/N 1+l/ as a func- 
tion of A^- 1 - ( /, 1 /(i-3^) for N = 64; 12 g ) 192, 256, and 

N = 512. In both cases, all data are described by a single 
master curve with \ 2 = 1-40 and 1.45, respectively. 

In view of the striking difference between planar and 
spherical geometry that we observe for the independence 
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FIG. 4: Average translocation time r (in units of Monte Carlo 
sweeps) normalized by N 1+v for chains escaping from (a) pla- 
nar confinement, as a function of the inverse driving force per 
monomer (A/x) -1 ~ (R/a) 1 ^ , and (b) spherical confinement, 
as a function of (A/i)" 1 ~ <j) 1 ^ 1 ~ 3l '\ The collapse of the data 
points for different chain lengths confirms the validity of Eqs. 
16al and J6bll . respectively. 



dence that follows the nonlinear lower bound established 
by Kantor and Kardar |l3| . The driving force affects the 
translocation time via a linear dependence on the inverse 
chemical potential gradient, as predicted on analytical 
grounds Thus, the average translocation time can 

be described by a simple scaling relation, r ~ N 1+ " /A/z, 
which results in a geometry- dependent power law when 
expressed in terms of the length scale of confinement 
distance between parallel plates for a planar geometry 
and cavity radius for a spherical geometry. 
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of the translocation time (Fig. [2J , it is remarkable that 
earlier work [{| found good agreement between numer- 
ical results for spherical confinement and a theoretical 
prediction based upon Eq. (J2J and Eq. Q, i.e., a driv- 
ing force obtained from the free energy of uniaxial (i.e., 
planar) confinement. We ascribe this finding to the for- 
tuitous cancellation of two errors. Indeed, in Ref. |9j the 
translocation time is predicted to scale as 



, (7\ V" 



si 



V(3") 



(9) 



Combination of the linear N dependence of Eq. 
and the inappropriate expression for the free energy of 
confinement yields an overall chain-length dependence 
JV 1+1 /( 3 ^, so that, at fixed concentration, r is predicted 
to scale as TV 1567 , which coincidentally is in approximate 
numerical agreement with the lower bound r ~ N 1+v 
[Eq. |J3J]. In fact, since the data in Ref. [t| exhibit a 
collapse when scaled by the N dependence of Eq. © 
we conclude that those data corroborate our findings. It 
is more difficult to reconcile our findings in Fig. [21 with 
the apparent confirmation of the concentration depen- 
dence of Eq. ©. However, we note that the evidence 
in Ref. consists of a scaling collapse rather than an 
explicit numerical analysis of the power law. When per- 
formed over a narrow range of densities, such a collapse 
can be relatively insensitive to the precise power of <$>. 

In summary, we have investigated the translocation of 
a flexible polymer chain through a hole, when the driv- 
ing force is generated by confinement of the polymer. To 
clarify the role of the confinement geometry, we have con- 
sidered the escape of a polymer from a planar as well as 
a spherical geometry. For both cases, we demonstrate 
that the translocation time has a chain-length depen- 
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